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ABSTRACT 

A double inflationary model with enhanced power at large scales — Broken Scale 
Invariant perturbations (BSI) — leads to a promising scenario for the formation 
of large-scale structure. We use high-resolution PM simulations with a model for 
the thermodynamic evolution of baryons. Galactic halos are identifled as peaks 
in the distribution of mass points which had time to cool providing a physical 
biasing mechanism. The clustering properties of galactic halos are in agreement 
with both small-scale and large-scale characteristics of the galaxy distribution, as 
variances of counts in cells, 'galaxy' correlation functions, 'cluster' abundances and 
'cluster-cluster' correlation functions, and velocity dispersions at small scales. 



1. Introduction 

During the last years, increasing evidence has been accumulated that the stan- 
dard CDM model in a = 1 universe as suggested by inflation leads to inconsi- 
stencies. First, the APM survey^ showed more power on scales larger than the galaxy 
correlation length. Later, the COBE measurements^ of the microwave background 
anisotropy provided evidence, that the CDM model would require an antibiasing for 
conformity with the observed galaxy clustering. Later, a large number of observati- 
ons came in difficulties with the CDM model. Therefore, additional admixtures of 
hot besides cold dark matter, or a cosmological term, or an open cosmological model 
have been proposed. We took up the old idea of a multi-inflationary scenario"^. The 
model decouples the structure formation on large and small scale, producing a power 
spectrum with a break at a characteristic scale"*. In addition to the amplitude, it 
requires two parameters: flrst, the break height being proportional to the ratio of the 
effective masses during the inflationary evolution; second, the epoch of the transition 
between the inflationary phases which is related to the initial energy density of the 
inflaton fleld. This epoch determines the scale length Is = 2tt jk^ of the break. 

We consider a special double inflation spectrum, which — using the linear theory 
— provides a good description to the large scale matter distribution in the universe^, 
also the APM angular correlation function can be well fltted^. A description of 
the primordial potential spectrum is given by the approximate formula (it does not 
describe the typical oscillations generated during the intermediate power law stage): 
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with A;^ = (27r/24)/i"^Mpc , A;/ = e^^h~^Mpc . We fold this spectrum with the 
standard CDM transfer functions (using a dimensionless Hubble constant h = 0.5) to 
get the power spectrum shown in Fig. 1. In the left we compare this BSI spectrum 
(sohd line) with the standard CDM model (dashed line), a ACDM (A = A/SH"^ = 0.8, 
= 0.2, dotted line) and a MDM model (Oqqj^ = 0.7, 0^ = 0.2, 0^ = 0.1, dash- 
dotted line). The standard CDM spectrum has more power at the scales relevant for 
galaxy formation, while BSI mostly resembles MDM at large scales, but has more 
power on small scales, almost as high as ACDM. 
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Figure 1: Power spectra (left) and variances of 'galaxy' count in cells. 



We have simulated the formation of large-scale structure employing a particle- 
mesh code'' and using a resolution of 512"^ cells and 256"^ particles for BSI®, compa- 
risons are made for COBF-normalised CDM with 256"^ cells and 128"^ particles. The 
box sizes lie between 25 /i~^Mpc and 500 /i~^Mpc to get a wide dynamical range (cp., 
Gottlober, this volume). On the right of Fig. 1 we compare simulated count-in-cell 
variances at 2; = as function of radius with observations in the IRAS catalogue^ 
(triangles) and the Stromlo-APM survey^° (circles). There the BSI curve (solid line) 
has been shifted vertically by a "biasing" factor 6 = 2 in order to normalize fluctua- 
tions to observations at 8/i~^Mpc , while the CDM curve (dashed line) has b = 0.9. 
On scales (15 — 50)/i~^Mpc, the CDM fluctuations are too low (outside of the 2-cr 
error range), at larger scales increase the errors. In contrast, the BSI simulations fit 
the data quite well, thereby testing the slope of the proposed spectrum. 

The used code includes a model for thermodynamic evolution of baryons (heating 
by shocks during the first anisotropic collapse, the standard cooling of the primordial 
matter, and a reheating of cooled particles to 5 X 10"* K with a fixed propability of 
0.85 of the baryons, say by supernovae, cp. Kates, this volume). The thermodynamic 
information is attached to the dark matter, not allowing a description of genuine 
hydrodynamic effects. We get an estimate of the local gas temperature used to 



describe the condensation of baryons which helps to identify different classes of objects 
in the universe. 

2. Galaxy halos 

After the first stage of pancake formation and growth, the particles virialize in 
potential wells possibly giving rise to galaxies. The galaxy finding procedure used 
is a modified density peak prescription. Based on previous experience'', we search 
for local maxima in the density contrast of cold particles exceeding a threshold 6th- 
We expect that shallower maxima do not lead to galaxies. The mass of 'galaxies' 
is determined by summing the contributions of all particles within a characteristic 
radius. For each simulation, several halo catalogues were constructed using thresholds 
6th = (1-5 - 3)0-5. 




Figure 2: 'Galaxy' mass spectra (left) and 'cluster' mass function (right). 



In Fig. 2 we show 'galaxy' mass functions obtained in a 25/i~^Mpc simulation at 
redshifts z = 2, 1. 5, 1,. 5,0 (curves from below). The dashed interpolation at 2; = 
corresponds to a Schechter fit 

$(M)dM = $,(M/M,)-*'exp(-M/M,dM, (2) 

with Af* = 3 X IO^^Mq , = 10~^ h^Mpc~^ and p = 1. These parameters are com- 
parable to the observations^^, but there, the mass spectra are steeper (the simulated 
CDM galaxy mass spectra of the same resolution are still much fiatter, providing too 
high numbers of heavy galaxies). This result is partly due to overmerging typical for 
our dissipationless DM simulations and the PM dynamics. However, the strong over- 
abundance of high mass halos is typical for the large power in the COBF-normalised 
CDM spectrum. We do not intend to make more detailed comparison with observed 
mass spectra of galaxies, but we use the identified halos with reasonable number den- 
sity to compare their clustering properties with the observed large scale structure of 
the universe. 



In the left of Fig. 3, we show the two-point correlation function for the galaxy 
catalogue constructed from a 200 h~^Mpc simulation (dash-dotted line). It is well 
described by a power law ^ = {r/ro)~^ for l/i~^Mpc < r < 15/i~^Mpc , with slope 
7 1.6 and correlation length Tq = 5.5/i~^Mpc (straight line), in agreement with 
the CfA-catalogue^^. The observed slope is influenced by peculiar velocity distortions 
as demonstrated by the dashed curve, extracted from all DM particles and imposing 
a linear bias b = 1.5. This is approximately the bias of a typical galaxy catalogue. 
The upper solid line is the correlation function of 'cluster' halos selected from the 
500/i~^Mpc simulation, which is enhanced by the same biasing factor; the slope and 
the correlation length Tq ~ 15/i~^Mpc are typical for cluster data (cp. Bahcall, and 
Bohringer, this volume). As in the case of the mass spectra, in the COBE-normalized 
CDM simulations the correlation function has a too steep slope (this is characteristic 
for low bias CDM^^^. 



10.0 
S 1 .0 
0.1 

1 10 1 10 100 

r [h^' Mpc] r [h^' Mpc] 

Figure 3: Correlation functions (left) and relativ-e peculiar velocities (right), 

3. Analysis of Galaxy Clusters 

Clusters of galaxies provide a sensitive test of cosmological scenarios. Their forma- 
tion is connected with the recent decoupling of large masses from the cosmic evolution. 
X-ray surveys^"^ and optical cluster catalogues are used to define a general 'cluster' 
mass function^"*. It tests the slope of the power spectrum at the break scale. Based 
on the Press-Schechter theory and using BSI primordial power spectra, we derived 
strong bounds on the break scale ks^^. To verify this conclusion we identify galaxy 
clusters in the 500/i~^Mpc simulations with a search radius of 1.5 /i~^Mpc ; it gives a 
reasonable cluster correlation radius and a positive correlation range up to 80 /i~^Mpc 
(cp. Fig. 3), like in observations. In Fig. 2 (right), we compare the cluster mass 
function with the data^"*, both the slope and the normalization fit the data (the solid 
line denotes BSI), only at Af < 6 X 10^^ Mq one notices a slight overproduction of 
cluster halos. Our estimate does not use any additional parameter besides the input 
perturbation spectrum, in particular it does not depend on the details of the thermo- 
dynamic assumptions in the simulations (we collect all particles in the cluster halos). 
The cluster formation is a quite recent process for the BSI spectra. In the contrary. 




the CDM model leads to a strong overproduction of lialos (cp. the dashed line), the 
mass function is to steep at cluster scales, and clusters are formed at high redshifts, 
contrary to observations. 

The velocity dispersion of 'galaxies' at different scales is an important test of the 
form of the power spectrum. We show the relative peculiar velocities as function of 
the distance on the right of Fig. 3. The dashed and dotted lines are the mean and 1-cr 
deviations of the velocity dispersion as determined by the anisotropy of the small scale 
velocity dispersion^^. The three simulation curves correspond to 'galaxy' catalogues 
with different thresholds corresponding to mass limits M > (5, 8, f f) IO^^Mq (from 
below). The different velocity dispersions are due to the higher probability of finding 
large mass halos in dense clusters, both due to merging and over-merging. As in the 
case of the correlation and mass functions, the studied CDM model leads to much 
higher small scale velocity dispersions®. The simulated velocity seggregation as shown 
in Fig. 3 represents a sensible test of the model. 

The primordial BSI power spectrum provides definitive predicions for the galaxy 
clustering. Our 'galaxy' selection procedure should be compared with full hydrodyna- 
mical simulations to test its trustworthiness. Certain statistics of large-scale structure 
as the correlation functions and the velocity fields are insensible to this scheme. They 
show that the BSI model passes many cosmological tests. The most important que- 
stion as compared with alternatives as CHDM or ACDM is the redshift dependence 
of clustering, and in particular, a quantitative comparison with predictions of the 
development of the halo mass function as discussed in Fig. 2. 
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